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Abstract 

We analyse the behavior of the density of states in a singlet s-wave superconductor with weak magnetic impurities 
in the clean limit by using the method of optimal fluctuation. We show that the density of states varies as In N(E) oc 
— \E — Ao| ( - 7 ~ d - >/ ' 4 near the mean field gap edge Ao in a d-dimensional superconductor. The optimal fluctuation in 
d > 1 is strongly anisotropic. We compare the density of states with that obtained in other recent approaches. 



Studies of spectral properties of disordered super- 
conductors remain an active area of research as they 
help advance our understanding of the competition 
between disorder and interactions. Recent years wit- 
nessed a renewed interest in the behavior of the density 
of states (DOS) in singlet s-wave superconductors with 
magnetic impurities. Such impurities are pairbreaking, 
and are characterized, in the weak scattering limit, by 
the spin flip scattering time, t s . In the self-consistent 
Born approximation (SCB) [1], the dimensionless pa- 
rameter controlling the suppression of the single parti- 
cle spectral gap, Ao, is Ar s , where A is the amplitude 
of the superconducting order parameter. In this paper 
we consider the clean limit, Ar s 3> 1, where the SCB 
approach yields a finite spectral gap, Ao ~ A, with the 
DOS N(E) = at energies E < A . 

It was argued in Ref.[2] that rare regions where lo- 
cal impurity concentration is high enough to locally 
destroy superconductivity lead to a finite density of 
states at the Fermi level. The argument was similar to 
the method of optimal fluctuation (OF), well known 
from studies of doped semiconductors [3] : in averaging 
over all the realizations of the impurity distribution, 
the probability of finding the realization which locally 
destroys the gap determines the residual DOS. Later, 
Lamacraft and Simons [4] considered in detail the en- 
ergy dependence of the DOS below the mean field gap 



in a dirty superconductor, where the scattering rate 
due to potential scattering greatly exceeds the spin-flip 
pairbreaking scattering rate, 1/t s . 

Very recently we analysed the subgap DOS in a clean 
s-wave superconductor, where the spin-flip scattering 
is dominant [5]; we argued that at least in some cases 
this limit is relevant experimentally. Here we briefly re- 
view the results of Ref . [5] and then present a more de- 
tailed comparison of the DOS obtained in Refs. [2,4,5]. 

We consider a mean field hamiltonian 



H = £t 3 + A(r)na 2 + U, 



(1) 
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where £ = -V 2 /(2m) — /i, fi is the chemical potential, 
Ti and o"i are the Pauli matrices in the particle-hole and 
the spin space respectively. The potential due to mag- 
netic impurities U = U(r) • s, where s is the electron 
spin operator, U(r) = JS;<5(r — n), J is the ex- 
change constant, and Si is the impurity spin at a site i. 

For an energy, E < Ao, OF is the most probable 
configuration of impurities that creates a state at E, 
and therefore contributes the most to the DOS [3]. OF 
provides nonperturbative corrections to the DOS de- 
termined in the framework of SCB. In essentially all the 
energy range below the gap the size of the OF is signif- 
icantly greater than the distance between impurities, 
so that the exact impurity potential can be replaced by 
a smooth function, and its probability density is well 
approximated by a Gaussian [3] with a width Uq = 
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riimpJ 2 S(S + l)/3, so that r s 1 = 2tyNoUq, where n imp 
is the impurity concentration, and Nq is the normal 
state DOS. The DOS is then given by In [N(E)/N Q ] w 
— S[U opt], where S[U op t] is obtained by minimizing 



S[U] 



2XJI 



d d rU 2 (r) + A I EfUl - E 



(2) 



with respect to both the potential U and the Lagrange 
multiplier A. We assume OF to be ferromagnetic [5], 
so that E[U] = <*|fr 3 + Aon + [7|*} where * is the 
normalized spinor wave function of the particle in the 
OF, and U(x) = -\US(^*{x)^{x)) is now scalar. 
In d = 1 we find (e = \E\/A ) 



U(x) 1-e 2 

2A e + cosh(2iVl - e 2 /Co) ' 
So = S[U opt ] = 87r(A r s ) \y/l-e 2 



e arccos e 



(3) 



(4) 



where fo — vp / Aq is the coherence length and vp is the 
Fermi velocity. The size of the OF, L x ~ fo/vl — e 2 > 
£o, and the depth of the potential \U\ < 2Ao. In partic- 
ular, near the gap edge, 1 — e <C 1, when [!7| <C Ao and 
L x > Co, we find S ~ (87r/3)(A r»)(l - e 2 ) 3/2 , which 
allows for a simple interpretation. In an OF of depth 
U and size L x the energy of the bound state is E ~ 
U + Ao + v F /L 2 Ao, subject to optimization of Eq.(2). 
This means \E — A | ~ vp/L x A ~ |{7|, from which 
an estimate for So ~ L x U 2 /Ui ~ (A t s )(1 - e 2 ) 3/2 
follows immediately. 

In d > 1 and for 1 - t < 1 we must compare the ac- 
tion for the isotropic OF with L ~ Co /V 1 — e 2 (when 
kinetic energy is vf/L), and for an anisotropic OF with 



L, and transverse size Lt 



a/a 



(corre- 



sponding to a wave function ^(r) = exp(ifci?£E)$(a;, y), 
where $ is a slowly varying function, so that 
l/(mL 2 ) ~ vp/L x ). The anisotropic OF is favored by 
a factor (E F /A ) (d - 1)/2 (1 - e)"^" 1 '/ 4 and gives [5] 



find S D /S iao ~ (A t s ) 2/3 (1 - e )< d ~ 2 )/ 4 , which, at the 
crossover, gives (A t s ) (10 " m)/ ' 6 > 1 for d < 3. Con- 
sequently, the OF discussed here still yields a higher 
DOS than that for purely diffusive motion. There- 
fore we expect that the structure of the OF near the 
crossover between the ballistic and diffusive regimes 
still resembles closely that given by us above, and that 
our results remain at least qualitatively valid there. 

Balatsky and Trugman [2] considered only the DOS 
at E = 0, where they needed a large volume fluctua- 
tion, V > £ d , which is less probable and yields lower 
DOS than that of Eq. (5). However, in the same spirit 
we should investigate whether local suppression of the 
gap from Ao to E due to a large number of impurities 
with uncorrelated spins (as opposed to a ferromagnetic 
OF above) is advantageous. For 1 - e < 1 the local 
pairbreaking rate, 7, needed to reduce the gap to E is 
7t s ~ 1 + (1 — e)(AoT 3 ) 2//3 , and the volume of the re- 
gion has to be at least equal to that of the anisotropic 
OF to avoid high kinetic energy cost (this is an under- 
estimate since it ignores proximity coupling to bulk). 
In that case we obtain the optimal action Sbt/So ~ 
(Aot s ) 1/ ' 3 (-Ef/Ao)c, where c = n imp X F is the atomic 
concentration of impurity atoms. As a result, for real- 
istic values of c and clean samples Sbt So, and the 
DOS given by the action in Eq.(5) is higher. 

We conclude that the action obtained in our ap- 
proach gives the most optimal fluctuation and the high- 
est DOS compared to other analyses. 
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We now compare this result with that of Ref.[4] 
for the diffusive propagation of the states within 
OF. It is clear that even in the clean case consid- 
ered here, with dilute magnetic impurities, a transi- 
tion to the diffusive regime occurs when the size of 



OF L > v F T s , or 1 



< (A 7 



In our nota- 



tions the result of Ref. [4] for Aot 3 3> 1 reads Sd = 
(Aor s ) 5/3 (£ F /A ) d " 1 (l - e) ((i ~ d)/i . Consequently, at 
the crossover point the action from Eq. (5) is smaller, 
S D /So ~ (Ep/Ao) {d - 1)/2 (A r s ) 7 ^ » 1, and the OF 
found here corresponds to a greater DOS. As the size 
of the OF increases even further, the anisotropic fluc- 
tuation becomes insupportable due to diffusive mo- 
tion. Nonetheless, even for an isotropic fluctuation, we 
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